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Abstract 

We explore further the Hamiltonian formulation of Yang-Mills theory in 2-1-1 dimensions 
in terms of gauge-invariant matrix variables. Coupling to scalar matter fields is discussed 
in terms of gauge-invariant fields. We analyze how the screening of adjoint (and other 
screenable) representations can arise in this formalism. A Schrodinger equation is then 
derived for the gluelump states which are the daughter states when an adjoint string breaks. 
A variational solution of this Schrodinger equation leads to an analytic estimate of the 
string-breaking energy which is within 8.8% of the latest lattice estimates. 



1 Introduction 



Yang-Mills gauge theories in two spatial dimensions, especially their nonperturbative prop- 
erties, are interesting for many reasons. First of all, they can be a model for the more 
realistic, but also more complicated, (3 + l)-dimensional theories. They have nontrivial 
dynamical content and propagating degrees of freedom in contrast to Yang-Mills theories 
in 1 -|- 1 dimensions, yet they are more amenable to mathematical analysis with some of 
the more recent techniques than their (3 + l)-dimensional counterparts. Beyond being a 
testing ground for nonperturbative techniques, they are also of direct relevance in, at least, 
one physical context, namely, Chromodynamics at high temperatures. (Recall that in an 
imaginary-time formalism, all Matsubara modes are suppressed at high temperatures, ex- 
cept for the zero modes of bosonic fields, leading to a dimensionally reduced theory [1].) 
Yang- Mills theories in three (or 2-1-1) dimensions are also of interest from the point of view 
of gauge-gravity duality, particularly in the context of the recent proposal for string and 
gravity duals of the maximally supersymmetric Yang-Mills theory in three dimensions [2]. 
Nonperturbative techniques for strong coupling analyses can also help elucidate some of the 
details of the gauge-gravity duality. 

A few years ago, a Hamiltonian analysis of pure Yang-Mills theory in 2 -|- 1 dimensions 
was developed based on a gauge-invariant matrix parametrization of the gauge potentials 
[21 d]. This led to the calculation of a wave functional for the vacuum state of the theory 
[5]. (For reviews, and other approaches and related papers, see [Q [71 [8].) The vacuum 
expectation value of the Wilson loop operator could then be evaluated using this wave 
functional and gave the result 



where R refers to the representation for the Wilson loop, and Ac is the area of the loop C. 
The string tension aR was given as 



where e is the coupling constant, and cr, ca denote the quadratic Casimir values for 
the representation R and for the adjoint representation, respectively. The area law ([T]) is 
consistent with confinement of charged representations. The calculated value of the string 
tension is in very good agreement with lattice estimates |9l [10] . More recently, motivated 



{Wr{C)) = exp[-(T^j^c] 





(2) 
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by our Hamiltonian analysis, there has been an attempt to estimate glueball masses 
We should note however that the vacuum wave function used in [11], which is somewhat 
conjectural, agrees with ours, derived in |5j, only in the high and low energy limitsll] An 
interesting variant of the Hamiltonian analysis is an 'anisotropic' formulation with two 
different couplings corresponding to the two spatial directions [8] . The theory is solvable in 
the extreme anisotropic limit; however, the approach to the isotropic limit is not entirely 
clear. In any case, from the results quoted earlier, it is clear that it is worth pursuing and 
elaborating on our Hamiltonian approach. 

The analysis was done in what may be characterized as a continuum strong coupling 
analysis. (This is different from the strong coupling analysis on the lattice.) For Yang-Mills 
theory in two spatial dimensions, has the dimensions of mass and the expansion parameter 
is /k or k/e^, where A; is a typical momentum scale. Modes of momenta much smaller than 

have to be treated nonperturbatively, and can be analyzed in a k/e^ expansion, while 

modes of momenta much larger than can be treated perturbatively in an /k expansion. 

Continuum strong coupling refers to an expansion treating modes of small momenta as 

the dominant contributions to any physical quantity. The nature of the expansion brings 

up many questions immediately. For example, the area law ([T|) and the string tension ([2]) 

are obtained for any representation, including representations of zero A^-ality (say, for the 

gauge group SU{N)). But these representations can be screened. Interpreting the Wilson 

loop as the propagation of two heavy external charges, we see that a gluon can bind to the 

external charge producing a color singlet, if the external charges are in representations of 

zero A^-ality. (The color singlet state produced by binding gluons to the external charge 

is generally referred to as a gluelump.) Thus, when the separation of the external charges 

becomes large, the energy in the string of gauge fields connecting these charges becomes 

large and it is energetically favorable to pair-produce gluons which bind to the external 

charges, neutralizing their color charge. There is no potential between the gluelumps, since 

they are color-neutral, and so we do not expect an area law for the Wilson loop for these 

representations. The process of the state of the two heavy charges connected by the string 

making a transition to the state of two gluelumps is called string-breaking |12[ [T5] . In our 

previous analysis, we did not see string-breaking. There is no inconsistency, since we were 

^ As we emphasize in the Appendix, in carrying out the regularization of the Hamiltonian and the "local" 
operators on which it acts, the order in which the regularization parameter of the Hamiltonian and the point- 
splitting for the operator are taken to zero is important in preserving important symmetries such as Lorentz 
invariance. We expect that this is one of the reasons for the discrepancies between the two wave functionals. 
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using the continuum strong coupling expansion, and string-breaking would not be seen to 
the order we calculated, but we need to understand how screening can arise as we improve 
on this approximation. 

Closely related to this is the question of whether the formula ([2]) for the string tension 
can be exact in any sense. Are there corrections to the formula even for representations of 
nonzero A^-ality? There are arguments in the literature, based on the 1/A^-expansion, that 
the ratios of string tensions for different representations should deviate from the ratios of 
the Casimir invariants [14J. The breaking of the string for a screenable representation is a 
transition from a single color-singlet operator to two color singlet operators. As such, the 
transition matrix element is suppressed at large N. So, one possibility is that the formula 
for the string tension, for the non-screenable representations, is exact in the large N limit. 
This question was numerically analyzed in detail recently by Teper and Bringoltz [lOj. They 
find that, while the formula ([2]) for string tension differs from the lattice result by only about 
0.88% as — > oo, the deviations are still statistically significant. It is therefore important 
to understand and calculate possible corrections. The approximation of the wave functional 
for modes of small momenta has a simple form, a Gaussian function in the appropriate 
variables; this simple form was used in evaluating the vacuum expectation value ([1]). While 
the closeness of the agreement between the lattice results and the values given by ([2]) shows 
that the Gaussian approximation is a very good first step, the wave functional has non- 
Gaussian terms and, therefore, it has the potential for explaining corrections to the formula 
for the string tension. This is an issue that deserves further analysis. 

In this paper, we will address the problem of screening of representations of zero A^-ality, 
deferring the issue of corrections to the string tension to a later publication. As mentioned 
earlier, screening is due to the formation of color singlet gluelump states, where gluons 
which are pair-created from the vacuum bind to the external charges. Here we will argue 
for the possibility of such a bound state from the higher non-Gaussian terms of the wave 
functional. We shall then construct the gluelump state, treating the external charge as 
a heavy scalar in the adjoint (or some other screenable) representation. A Schrodinger- 
type eigenvalue equation will be obtained, within certain approximations explained in the 
text. We then carry out a variational estimate of the gluelump energy. This gives us the 
estimate for the string-breaking energy which we compare to the value obtained from lattice 
simulations [TS]. 
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For carrying out this analysis, we need to develop a gauge-invariant Hamiltonian ap- 
proach for scalar fields coupled to the Yang-Mills fields. This will be the subject of the 
next section. This formalism is also the first step in extending the continuum strong cou- 
pling analysis of the gauge theory to the case with matter fields added, and as such, it 
is of interest beyond understanding the physics of screening. For instance, it can lead to 
an analysis of the supersymmetric theory, once we have further extended this to include 
fermion fields. We may note that there has been a recent proposal about the gravity dual 
for (2-|-l)-dimensional Yang-Mills theory with sixteen supercharges [2]. The gauge theory in 
question can be regarded as the dimensional reduction of AA = 4 supersymmetric Yang-Mills 
theory from R x S"^ to R x S^. A prediction for the leading strong coupling contribution 
to the masses of the operators built out of the scalars fields has also been computed from 
the string theory side [2J. Since the gauge- invariant framework that we shall develop in the 
paper naturally leads to a strong coupling expansion for the gauge theory in the continuum, 
there is the potential for obtaining similar results directly in the strong coupling regime of 
the gauge theory and for the possible confirmation of some details of the gauge-gravity 
duality. 

We also note that three-dimensional Yang-Mills theory with matter fields in the adjoint 
representation naturally arises at finite temperatures in an effective theory of deconfinement 
|15j . The "matter fields" in this the Polyakov loops around the imaginary time 

direction. They take values in the gauge group; nevertheless some of our analysis can be 
adapted to analyze that theory. We consider this paper again as a prelude to such an 
analysis. 

In section 2, we construct the Hamiltonian for the (2-|- l)-dimensional Yang-Mills theory 
coupled to adjoint scalars in terms of gauge- invariant variables. In section 3, we give 
qualitative arguments for how screening could arise in the Hamiltonian approach. In section 
4, we then use the gauge-invariant framework outlined in section 2 to work out the action 
of the Hamiltonian for gluelump states of very short spatial extension, appropriate to an 
extreme strong coupling approximation. This gives a lower bound on the string breaking 
energy. The extension of this to gluelump states of finite spatial extent is given in section 
5 where a Schrodinger-type eigenvalue equation is obtained. In section 6, based on this 
eigenvalue equation, we obtain an estimate of the string-breaking energy and compare it 
with available lattice data. There is a short summary/discussion section and the paper 
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concludes with an appendix with some comments on how regularization has to be carried 
out for various terms in the Hamiltonian. 



2 Scalars coupled to the YM field 

In this section, we shall set up the essentials of coupling a real scalar field to Yang-Mills 
theory in two spatial dimensions. We will consider an SU (N) gauge theory and the scalar 
field will be in the adjoint representation of SU{N), although various formulae can be easily 
generalized to any representation. 



The action for the theory is 

S = — \r [ d?x Tr(F2) + - I d?x Tr {DQct)f - {Dcjyf - M'^cf^ 
2e^ / 2 



(3) 



The corresponding Hamiltonian is given by 

n = 



(fx 



2e2 
6^ 



+ 



+ 



2e2 



+ 



+ 2{Dcj)y{D<l)r 



(4) 



In writing this expression, we have made the gauge choice Aq = 0. A = ^{Ai + iA2), A = 
^{Ai — iA2) are the two complex spacial components of the gauge potential. D, D are the 
holomorphic, antiholomorphic covariant derivatives 



D(^ = d(j)+ [A, 



(5) 



In our conventions. 



A = -iAH" 



where the SU{N) generators are normalized such that 



(6) 



(7) 



We would like now to rewrite the Hamiltonian in terms of gauge-invariant variables 
following the approach used in O SI E] for pure Yang-Mills theory. As before, we use the 
parametrization 

A = -dMM~^, A = M'^~^dM^ (8) 
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If the gauge transformations take values in the Lie group G, M is a complex matrix taking 
values in G^, the complexification of G. For the case of G = SU{N), which is what we 
shall consider in this paper, M G SU{N)^ = SL{N,C)- 

Time-independent gauge transformations are realized as left rotations on the matrix M, 

M{x) g{x)M{x) , g{x) G SU{N) (9) 
As a result, the theory can be entirely rewritten in terms of the gauge-invariant variables 

H = M^M , x = M^(pM^~^ (10) 



Before we derive the expression of the full Hamiltonian ^ in terms of these gauge- 
invariant variables, it is useful to first review the case of the pure Yang- Mills theory. The 
simplification of the Hamiltonian in terms of the variables M, is as follows. Introduce 
the right-translation operator pa for M and the left-translation operator pa for by 



[pa{x),M{y)] = M{y){-ita)5^^\x-y) 
\pa{x),M\y)] = {-it,)M\y)5^^\x - y) 



(11) 



It was shown in [3] that in terms of these, the kinetic energy operator T may be written as 

.2 r ^2 „2 



T 



e I 0" _ e' 

Qar {x, u)Kab{x)Qhs {x, v) 



Urs{u,v)Pr{u)Ps{v) 



(12) 



where Kab = 2Ti{taHt},H ^) is the adjoint representative of and 

1 



Gma{x,y) 
Qma{x,y) 



tt{x - y) 

1 



6ma-e'\^-y^^/^{K{x,y)K~\y,y)), 

5ma-e-\^~y\'''{K-\y,x)K{y,y))^ 



■k{x - y) 

These are the regularized versions of the corresponding Green's functions 

1 1 



G{x,y) 



Gix, y) 



(13) 



(14) 



7r(x - y) ' ' 7r(x - y) 

The parameter controlling the regularization, e, acts as a short-distance cut-off. Expression 
()12|) is to be used on functionals where the point-separation of various factors is much larger 
than ^/e. (See Appendix for more details on regularization.) 
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As discussed in [3] there is an ambiguity in the parametrization ([8]) : (M, M^) and 
{MV{z), y(z)M"l") give the same gauge potentials A, A, where V, V are, respectively, 
antiholomorphic and holomorphic in the complex coordinates z = xi +1x2 and z = xi —1x2- 
As a result, all physical observables in the theory should satisfy the so-called holomorphic 



mvariance 



H V{z)HV{z) (15) 
The particular choice of regularization used in (jl3p respects this holomorphic invariance. 



When the theory is rewritten in terms of gauge-invariant variables, Pa,Pa can be realized 
as translation operators on H = M, 



[paix),Hm = H{y){-ita)5^^\x-y) 
[Pa{x),H{y)] = {-ita)H{y)5^^\x-y) 



(16) 



In OH], we have argued that the wave functions can be taken to be functionals of 
the current J = {cA/Tr)dHH~^ , where ca is the quadratic Casimir invariant defined by 
0^6""^ = xhe action of the operators p, p on the current J is expressed in terms 

of the commutation relations 



[Psiv], Ja{z)] 
\pr{u), Jb{w)] 



.CA 



i — Kas{z)dz5{z - v) 



vr 



(17) 



where V^j ah = ^d^Sab + ifabcJciw). Using (fT2|) and (fTTl) . the action of T on wavefunctions 
of the form ^{J) can be expressed as 



Tym ^{J) 



27r 



a;"(z)^-7- 



+ 



6 



6J^(z)6J^{w) 



^(J) (18) 



where 



Ara{w,z) = -[dzllrs{w,z)] KJ{z) 
^ab{z,w) = 'D^brKa(w,z) 



(19) 



The above expressions have been computed in [3] where we found that for small e, T 
can be further simplified as 



Tym^{J) = m 



Ja{z) 



SJa{z) 



+ j {V^G{z,w)) 



ab 



5Ja{w) 6Jb{z) 



^(J) + 0(e) (20) 



8 



where m = c^ca/Stt. 

As for the potential energy term, it is trivially checked that we can write 

Vym = I cPx = [ Sx : BrBr : (21) 

J 2e^ mcA J 

The regularized form of this expression is 

Vym = ^[[ a{x,y-,X)dJaix){Kix,y)K'\y,y))abdJi,{y)-^^^^] (22) 

fTT'CA Jx,y ^ 

<y{x,y;X) = — 

ttA 

where u{x, y; A) is a regularized 5-function, A is the parameter of regularization, and we 
should take the limit where \/A ^ 1/e^. The operator U'^^{x,y) = [K{x,y)K~^{y,y)]'^^ in 
(|23p is a "holomorphic" Wilson line connecting the two 5J's at different points. U°'^{x, y) is 
such that the regularized expression for Vym satisfies holomorphic invariance. In considering 
the application of T on this, and on other expressions in the rest of this paper, the correct 
procedure is to preserve A » e as we take these quantities to zero. This is what was used 
in [HIS]; it is also briefly explained in the appendix. 

Since T as defined in (jl2p is valid only for separations much larger than -^e, we need to 
use A ^ e to define the potential energy term. Correspondingly, we can define the kinetic 
term also with a scale A by the equation [1] 

r(A) = r + ylog(26/A) Q 

Q = e j a{u,v;e)Krs{u,v) (^pr{u) - idJriu)^ Ps{v) (23) 

The understanding here is that I ^ e; this operator is to be taken to act on functionals of 
J with a point-separation much larger than y/e, with e ^ and \/I ^ 1 /e^ . 

We now turn to the case with the scalar field added. Evidently the term J SASA does 
not change and we can still use the expression (jl2p . The full kinetic energy operator T is 
now 

1 /• 2 1 f 6^ 



2 ./ ^A^bA^ 2 

o2 



^ I ^rs{n,v)Pr{u)Vs{v) -^J (24) 
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where x is the gauge-invariant field defined in (jlOp . The action of p, p on x is given by 



\pa{x),xbiy)] 
[Pa{x),xb{y)] 



-fabc Xc{y) 6'^^\x-y) 



(25) 



The gauge-invariant wave function, ^{J, x), is now a functional of both gauge-invariant 
variables J and x- In terms of its action on such functionals T can be written as 

6 6 1 r 6^ 



T = TYM + im [ A^w, ^r'^X^iw) 
where Tym and Ac^(tiJ, i*) are given in p8]) . ([19]). 



6x''iw) 5J<^[z) 2 J Sx^'^x" 



(26) 



Similarly the full potential energy term can be expressed in terms of gauge-invariant 
fields as 

= Vym + j [^f^dx^'iVxT + ^^^^^ 



(27) 



For purposes of easy reference, we shall now collect various formulae and write the final 
form of the Hamiltonian for the Yang-Mills field coupled to adjoint scalars as 



r _ 1 

n = — j Ilrs{u,v)pr{u)ps{v) - - 



+ Vym + 



-dx\vxr-r^^ 

CA 2 



m 



5 



6 



+im I A,rf(u;,z)rV(^) 



5 6 



+ 



mcA 



: dJ^dJ" : 



6xHw) 5J<^{z) 

M-^X^X" 



(28) 



2 J Sx'^Sx" J \CA 2 
It is evident from our derivation that this formula for the Hamiltonian is trivially general- 
izable to the case of several scalar fields. 



3 Screening: a preliminary qualitative analysis 

The next logical step is the application of the gauge-invariant formalism of the previous 
section to the gluelump state. But before doing so, it is useful to go over some qualitative 
considerations of how screening of the adjoint and other representations of zero A^-ality could 
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arise in our Hamiltonian approach. As discussed in [5], the computation of the expectation 
value of the Wilson line in pure Yang- Mills theory is of the form 

{Wr{C)) = I dfi{A/g.) e-^(^) Wn{C) 

= J d^{H)e^'^^--^^'> e-^(^) Wr{C) (29) 

where dfi{A/g^,) = dfi{H)e'^^'^^'"'"'^^^ is the gauge-invariant volume element for Yang-Mills 
fields and e~^^^^ is given in terms of the vacuum wave function as e~^^^^ = ^'q^'o- Su)zw{H) 
is the Wess-Zumino-Witten action for the hermitian field H and d^{H) is the Haar measure 
for H viewed as an element of SL{N,C)/ SU{N). Explicitly, 

Sn^UH) = Tv{dHdH-^) + j e^^^^TiiH-^d^HH-^d^HR-^d^H) (30) 

The expression for the wave function shows that S{H) has the form 



47r2 

S{H) = -^-^ / dJa 



dJa 



{m + y/m? — V^) 

-Vabcj f^''\x,y,z)Ja{x)JM-W) + 0{J^) (31) 

The function f^^\x,y,z) has been given in [3]. 



Equation ([29]) shows that expectation values can be considered as averages in a two- 
dimensional field theory with the action S{H). In evaluating the string tension and obtain- 
ing the formula in terms of the Casimir values of the representations, we have basically used 
the low-momentum limit of the quadratic terms (terms which are quadratic in the current 
J) in this functional. These terms correspond to 

SiH) = ^ld'x F^F?^ (32) 

where = me^ = e^CA/'^.'K. In this approximation, the expectation value reduces to 
the average calculated using a Euclidean two-dimensional Yang-Mills theory with coupling 
constant g. There are many ways to evaluate these averages, one simple way is to go back 
to the original variables, the A's, and calculate, for example, in a gauge Ai = 0. In this 
gauge, the commutator term in Fij is absent and we have a Gaussian functional integral. 
This leads directly to the formula for the string tension. 

As mentioned in the introduction, the values we find are in very good agreement with 
the lattice analysis [9l [10], although the lattice results are accurate enough to show a 
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statistically significant deviation, even at large A'^, where the agreement is the best. Also, 
it has been argued, on the basis of the 1/A^-expansions, that deviations from the Casimir 
scaling for string tensions of different representations are possible |14] . Further, we know, 
on general theoretical grounds, that the Wilson loops in representations of zero A^-ality can 
be screened and will not show the area law behavior for large enough loops. Therefore, 
we turn to considerations of possible corrections to this result. The corrections can arise 
from the higher terms involving three or more powers of the current J. The string tension 
is defined by the limit of large Wilson loops. Therefore we may think of "integrating out" 
the high momentum modes in the two-dimensional theory defined by S{H) to obtain an 
effective low- momentum action Seff{H) and then using this to calculate the average of the 
Wilson loop. It is then easy to see that there are two types of corrections possible. The 
first set would correspond to a corrected quadratic term, due to terms which modify the 
"propagator" for the currents. This is equivalent to a modified coupling g^, in the effective 
two-dimensional Yang-Mills theory. Since this is at the level of the action and we have not 
yet introduced the Wilson line, there is no dependence on the representation of the Wilson 
line; it depends only on and numerical factors. Therefore, these corrections can affect 
the value of an, but the ratio aR/ap is unaffected. 

The second set of corrections involve the and higher vertices with the currents con- 
nected to the Wilson line via propagators. These cannot be reduced to propagator correc- 
tions and hence we can get representation-dependent modifications. In particular, screen- 
ing effects have to come from such terms. In principle, all the and higher vertices can 
contribute. To elucidate the nature of these terms, consider the contributions from the 
cubic vertex [5j. Written in terms of the original variables, the cubic term is of the form 
Th:{D^'^B[B,DB\) (Since we have phrased the calculation of the string tension in terms of 
the ^'s, we use the same variables for this argument.) Expanding in powers of A, we 
see that 

D-^ = d-^ -d-^Ad-^ + ■■■ (33) 
In the Ai = {) gauge, the cubic term is thus of the form 

Sfat - [(^~' - d'^Ad-^ + ■■■) diA2] " diA\ didA^ (34) 

Because of the overall Euclidean invariance, the Wilson line can be visualized as the propa- 
gation of two heavy charges. Consider then the insertion of the interaction term (j34p . The 
contractions of the A's from the expansion of will produce a sequence of terms which is 
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s 














-J 





+ • 




A set of terms whose sum can be viewed as the propagation kernel of a bound state 

(shown as thick gray hne) 



effectively like the iteration of the kernel in a Bethe-Salpeter equation for a bound state. In 
other words, we get terms which can be interpreted as the propagation kernel of a possible 
bound state of the gluons with the external charge of the Wilson line. For the nonscreenable 
representations, it is not meaningful to think of the bound state in isolation since it is not 
a color singlet. But for the screenable representations, a color singlet bound state can be 
formed. The picture which emerges is one where gluons are created from the vacuum and 
then bind to the external charges to form color singlets. This can lead to breaking of the 
single string connecting the external charges into two such singlet bound states. 

As for the quantitative calculation of these corrections, we know from the lattice data 
that they are small, except for the screening effects. It is, of course, important to calculate 
the corrections and provide an analytic justification for why they are small. The calculation 
of these corrections and the analytic justification for why they are small will be reported 
in another paper. Here we will pursue the analysis of string breaking for screenable rep- 
resentations. There are two related issues for the analysis of this question. There is a 
transition matrix element between the state of the single string connecting the external 
charges and the two singlet bound states. This matrix element will describe the process of 
the string-breaking, and, from general arguments, it should be suppressed by powers of N, 
in a large expansion. The calculation of this matrix element is quite involved and we will 
not carry this out here. However, since the energy at which the string breaks corresponds to 
twice the energy of the singlet bound state, we can calculate this string-breaking energy by 
calculating the energy of the singlet bound state of a gluon with the external heavy charge 
(whose propagation is one side of the Wilson line). This bound state has generally been 
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referred to as the gluelump. In our language, it corresponds to a wave function of the form 

*G = I d^xd^y f{x,y) dr{x)[K{x,y)K~\y,y)r'x\y) ^0 

= Og ^0 (35) 

Here J represents the gauge-invariant gluon which has a mass m. In reference j^, we have 
discussed a constituent picture for the formation of bound states, such as gluebaUs, in terms 
of elementary constituents of mass m corresponding to J. The gluelump, as described by 
([35]) . is consistent with that picture. The function f{x,y) is then the wave function of the 
constituents in the bound state. 

We can now move onto a more specific and quantitative estimate of the bound state 
energy rather than the general qualitative arguments given above. For this, the strategy will 
be to act with the Hamiltonian on the state (j35p and obtain a two-particle Schrodinger-like 
equation for the function f{x, y) characterizing the bound state. The details of this analysis 
are given in the next section. 



4 The action of the Hamiltonian on ^ 



G 



We will first consider the (continuum) strong coupling limit where we use the state ([35 
with the vacuum wave function ^'q = ex.p{—^S{H)) approximated by 



exp 



exp 



2m 2 



1 



M 



"8^ J d'xF'- - I d'x 



(36) 



We have used the perturbative vacuum for the X'^^ld. In the large M limit, because 
^> BV, corrections due to the interaction term in BV are suppressed by powers of 1/M. 

With the approximation of ^'o as in (I36p . the action of the Hamiltonian on ^g can be 
written aj^ 

n^G = M^G + Og Hym^o + ^o(TOg) 
e2 f 

/ ^rs{u,v)[priu),VYM][Psiv),OG]^0 



-— I Urs{u,v)[Psiv),VYM]\Pr{u),OG]'i'0 



(37) 



^We take the Hamiltonian to be normally ordered for the scalar fields so that there is no zero-point energy 
contribution in (1371). 
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Since ^'o is the vacuum for the Yang- Mills part, TCym^o = 0. For the third term and 
the last term on the right hand side, it should be kept in mind that pr has nontrivial action 
on X as given in (p5|) . 

In evaluating (TOg), it is useful to compare it with the action of the Yang- Mills kinetic 
term, which we have evaluated elsewhere 0]. Consider 



Notice that Og is obtained from this by replacing the last dJ, namely dJ^{y), by x^iv)- 
The action of T on the state (|38p can be written as the sum of four terms as given below. 

T ^2 = 1+ 11+ IV 



where U"'^{x,y) = [K{x,y)K~^{y,y)]°'^. We will consider the evaluation of these terms for 
small separations |x — y| <C with |x — y| ^ ^/e. (This means, of course, that the the 
width of / which controls the average value of |x — y| is small compared to 1/e^.) Term 
/ then leads to a mass m for dJ and the Coulomb potential V(x, y). Terms // and III 
contain a normal-ordering term, with the p, p acting on the 9 J's. The other terms arising 
from these two are negligible for very small separations; these involve at least one power of 
\x — y\ and three or more J's. Finally the last term gives m for the final dJ{y). Thus we 




(38) 




(39) 



find 




[2m + V(x, y)] f{x, y) : dJ^{x)U^\x, y)dj\y) : + 0{\x - y\) 



(40) 



The Coulomb potential in (I40p is of the form 



V{x, y) = m Ein(|x — y|^/2e) 



m \og{\x — yp/2e) 



(41) 



where 




(42) 
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In comparing this with TOg, we see that we will get terms similar to / and II, but the 
terms analogous to /// and IV will be absent when dJ{y) is replaced by xiv) because of 
(|25p . In particular we find 

„2 



I' 



- / Ursiu,v)[priu),[Ps{v),dr{x)U'^''ix,m x'iv) f{S,y) 



m + V{x,y)]fix,y) dr{x)U-\x,y)x'iy) + 0{\x-y\) (43) 



The normal ordering contribution from the term similar to // is actually zero, due to 
color contractions. 

„2 



Urs{u,v] 



[Psiv),dr{x)]u''\s, y) [Pr{n),x\y)] 

+dr{x)[psiv),U''\x,y)] [p,(n),x'(y)]] fix,y) (44) 
Using the commutation relations (|17p we find that the first term in //' can be written as 



= - J Urs{u,v)\psiv),drix)]U'''{x,y)[priu),x\m fi^,^ 

= mf d,Araiw,z}[Kiz,w)K'Hw,w)]^j''\'iw) fiz,w) (45) 



Explicit calculation of A^^ in [4j shows that (see eq.(A3) in Appendix 



K{w,w)K-^{z,w)+ J2 {z-w)^{z-w)'^Fab{w) + 0{e) (46) 

A,B>1 



2tT€ 

where a = \z — tfJp/e. Inserting (j46p into ()45p we find that only the first term in (j46p might 
contribute a finite e*^-order term, but this is actually zero due to color contraction. 



(a) 



m 







-a/2 



2e7r 



[K{w)K-\z,w)]^^ [K{z,w)K-\w)]^j'^'x'H f{z,w) 



(47) 



The second term in (j44p can be similarly evaluated 

„2 



^ I Urs{u,v)dr{x)[ps{v),U-\x,y)] [Pr{u),x\y)] f{x,y) 



dr{x) [Ur.M X, y) - iirs{y; y)] r'f'^'K-'ix, m'^^iy, y)x'{y)f{^, y) 

(48) 

In writing (j48p we used the fact that 

[pM, u-\^. y)] = r'^'K'^'ix, m'^'iy, y) H'^ - ^) - K'^ - y)] 5{v - y) (49) 
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Explicit calculation shows that 

Ilrs{y;x,y) - ]lrs{y;y) = — > ^ 



vr ^ — ' nn\ 

n=l 



(50) 



Inserting (jSOp into (|48p shows that the contribution of -^-^^j) is, indeed, negligible for very 
small separations. So we find that 

(TOg)=/ [m + V(f,y)]/(x,y) aj"(x)[7"''(x,y)x'(y) + (51) 

Jx,y 



(52) 



Finally, the remaining terms in (I37p can be simplified using 

47r /■ - (5 

[r,yyAf] = 2m Wm + — / T^dJ— 

CA J OJ 



The term ImVyM is due to both p and ]5 acting on Vym] this is not needed for the terms 
under consideration. We then find 

„2 



4m 



Urs{u,v) {\Pr{u),VYM]\Ps{v),OG] + [ps{v) ,Vy M][Pr{u) , Og]} 

fix,y) ■.dr{x)U-''{x,y)x\y): + 0{\x - y\) (53) 



/ 


"_v|- 


/x,y 


2m 



The e — > limit of (|20p . which has been used to simplify the right hand side of (|52p. is 
adequate for the calculation of the right hand side of (|53p . Notice that the resulting factor of 
— V^/2m, along with the mass m in (I5ip . will reproduce the first two terms of the relativistic 
combination Vm^ — V^. The remainder is of order \x — y\. 

Collecting various terms and results, we get 



x,y 



M + m- — + V(x, y) 
2m 



/(x,y)| drix)U-\x,y)x\y) 

+ 0{\x-y\) (54) 



So far, we have not considered the addition of Q as in (|23p . As mentioned earlier, the 
Coulomb potential in (|5ip is of the form Ein(|x — yp/2e) w log(|x — yp/2e). From the 
definition of Q, 

Q J dr{xW\x,y)x\y) f{x,y) 

= 2A [ dr{xW\x,y)x\y) fix,y) (55) 
vr 7 
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where, in the first step, we have indicated the term which gives a nonzero result among 
the various terms generated by the action of Q. We have evaluated the action of Q on 
dJ°'{x)U"'^{x, y)dJ^{y) before [Ij; the nonzero term has the same structure as in (jSSp . except 
for the replacement of by dJ^ . Thus, the result in (|55p is exactly of the same form as 
obtained for dJ"-{x)U"-'^{x,y)dJ^{y). 

The addition of (e^/2) log(2e/A)Q to the definition of T thus leads to the result 



x,y 



V2 

M + m- + V{x,y,\) 
2m 



/(x,y)| dr{x)U^\x,y)x\y) *o 
+ 0{\x-y\) (56) 



where V{x,y,X) = m (Ein(|x — y|^/2e) + log(2e/A)) w mlog(|x — y|^/A). (Notice that if 
f{x,y) is of the form a{x,y,X), the Coulomb potential vanishes.) Thus, within the ap- 
proximations used for this continuum strong coupling expansion, we see that the state 
becomes an eigenstate of the Hamiltonian 7i with eigenvalue E if f{x,y) obeys the 
eigenvalue equation 



2 



M + m- — +V{x,y,\) 
2m 



fix,y) = E fix,y) (57) 



This equation shows that there is a constituent picture for the gluelump. The gauge- 
invariant gluon of mass m binds to the external charge via the Coulomb potential. Since we 
have arrived at this equation keeping only the F^-teim in (Yang- Mills part of) log ^'o, this is 
appropriate for spatial momenta small compared to m. The kinetic energy for the "gluon" 
is thus the nonrelativistic expression. As mentioned before, this is the beginning of a series 
which sums up to V — [5j. Further, we find only the Coulomb potential. The key 
question is whether these approximations are adequate. If we are interested in a prediction 
in the continuum strong coupling limit, it is adequate to keep the nonrelativistic kinetic 
energy term. We cannot say, a priori, whether this is sufficient for comparison with lattice 
simulations, since a direct correspondence of the kinematic regimes in the two approaches 
is not obtained. However, experience with the comparison of the string tension and glueball 
masses suggests that this is a good starting point. 

The question of the potential energy V(x, y) is another matter. From the action of Tym 
on dJ{x)U{x,y)x{y), we get the Coulomb limit of the potential. When effects due to the 
vacuum wave function are included, this potential must tend to the linear potential at large 
distances. We know this is so, from the existence of a nonzero string tension. However, 
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a direct calculation is difficult and we shall present a somewhat indirect argument in the 
next section. We will also use a slightly different definition of which is more convenient. 
On the lattice side, the potential energy is seen to rise linearly with distance, with a value 
of the slope or string tension which is in agreement with our calculations, up to a point, 
and then it flattens out, indicative of the breaking of the string. Thus, just at the point 
where the string breaks, the constituents of the singlet state (the external charge and the 
"gluon" ) are separated far enough to be in the regime of the linear potential. For a sensible 
comparison with lattice data, we must understand and incorporate the linear potential. We 
take up this task in the next section. 

However, already at this stage, it is worth noting that, if we consider the strong binding 
limit (when the separation \x — y\ is negligible) and also the strong coupling limit where 
— <^ m^, the energy in (|57|) can be approximated by M + m. The resulting value 
2{Eq — M) = 2m may be taken as a lower bound on the string-breaking energy. 



5 The argument for the Unear potential 

The strategy we will consider for extracting the linear potential is as follows. We consider 
a state of the form 

\d) = [ f{x,y) tl;^-{x)W-\x,y)<P^''{y) \a) (58) 

Jx,y 

where W{x, y) is the open Wilson line connecting points x and y, W{x, y) = V exp ^— 

and ip and cj) are heavy fields. \a) specifies the state of the gauge field and can be expressed 

as |a) = Oa\0), where Oa is constructed from the gauge fields only. 

The matrix element of e"^"^ for states of the form (j58p can be related to the expectation 
value of the Wilson loop; the latter can be expressed in terms of the area law and the linear 
potential extracted from it. We use two dissimilar fields tp and (p to avoid the possibility of 
mutual annihilation during the (Euclidean) time-evolution. 

The action for the theory is taken as 

S = j d^x[{D^)%DiPf + {D4>f{D(Pf-M^{ij''i,'' + r(t>l] + Sym (59) 
The Hamiltonian has the form 

H = HyM + 'Hscalar + "Hint (60) 
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where the last term involves interactions between ip and A and between (p and A. This 
Hamiltonian commutes with N^, the number of V-particles and A'^^, the number of (j)- 
particles, separately. Thus matrix elements of the Hamiltonian are block-diagonal and 
we can restrict attention to the = = 1 block consistently. (A real adjoint scalar 
particle can, in principle, annihilate with another such particle. This would be an added 
complication which is irrelevant for our approach to the linear potential. So we avoid it by 
choosing two distinct complex scalars to construct the state. This is not the minimal way 
of doing this. A minimal way would involve one complex scalar, but the final results are 
the same, so we use the more convenient choice of two complex scalars.) 



First consider the normalization of states of the form (j58p . For this, the states are at 
equal time. We can use iI){x)iIj^ {y)\a) = 6{x — y)\a), since Oa has no ip^s in it. We then find 



{a\a) = f*{x,y)f{x,y) {a\a) (61) 

Jx,y 

The evaluation of {a\a) can also be simplified if we approximate the vacuum wave function 
by the expression (|36p . We then find 



{a\a) = J dfi{A/G*) exp J d^xF"^ 



OaOa (62) 



We can choose the gauge Ai = 0, whereupon the action is quadratic in A2 with a propagator 
(Ai(x)^^(y)) = <^»''<5(xi - yi) J ^e'M-^'y^) ^ (53) 

The actual evaluation of (j62p can then be done by Wick contractions using (j63p . 

We now turn to the matrix element (/3|e~^'^|a), where is another state similar to 
Id), with Ojs and a function h{x,y) replacing Oa and f{x,y), respectively. 

(/3|e-«^|Q) = [ h*ix',y')fix,y) {my')Wiy' , x')i^ix') e"^- ijHx)Wix,y)^\y)\a) 

h* (f y')/(x, y) (0| [Op<t>{^)W{^, x')V'(x')] {x)W{x, y)cP^ {y)Oa \ 0) 

(64) 

The set of fields at the left end have been shifted by r in Euclidean time. Since the fields are 
heavy, we can write, for large M, ^(t') = e^"^ cjye^'^'^ ~ e~*^'^(/)(0), so that (i?5>(r, x')0^(O, x)) ~ 
can then simplify the matrix element as 

'Ht\~\ _ I r*/-j „-t\ffz: „-t\ „-2Mt 



|e"«n«) = / h*{x,y)f{x,y) e-''''^ (0| 0^(r) W{C) O,(0)|0) (65) 
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where 0^(0) = O^ir = 0) and W{C) is the Wilson loop W{C) = TrV exp{- §^ A). The 
path C is a rectangle connecting the points (x, 0), (y, 0), (y, r) and {x, r). The evaluation of 
the remaining expectation value in (j65p can be understood using Wick contractions. There 
are two types of contractions involved: 

1) Wick contractions within W{C) and separately within O/^Oa 

2) Wick contractions between W{C) and OjsOa- 

We neglect the second type of contractions; we comment on this later. In this approximation, 
we have 

(o| [Ofs]^ w{c) o^mo) « (0|O;3(t) o,(o)|o) {o\w{cm 

= e~^\'^~'^^ (0|O/3(r) O«(0)|0) 

= e-'^l^-s^l^ {(5\ e-^^^^^|a) (66) 

Once we approximate as in the first line of this equation, the expectation value of W{C) 
can be evaluated, by Euclidean invariance, as an equal-time average. Also, notice that only 
q-Hymt appears in the last line. Thus, within the approximations we have made. 



h*(l,W^ e~^" j f^^W(l)'^\a) = I h*f e-2^^"-"l^-^^"(/3|e-^-'^n"> (67) 
We may rewrite part of this integral as 

h*{x,y)f{x,y)e-''^''-y^^ = (0| [ / h*(PWi^]^^^ e-l--j7k[ / f^p^WcP^^JO) (68) 



ix,y 

Using this result, we finally get 



h*^Wi; e-^^ / f^P^W(|)'^\a) = {i]\e-(^M+a\x-y\)r / f^p^W(j)^e-'^'""^\a) (69) 



Since \P) is completely arbitrary within the A''^ = = 1 block, we conclude, by taking 
terms linear in r, 

n j fij^Wcp^a) = {2M + a\x-y\) j f'ilj^W(l)^a) + j fi;^W^^ TiYMOalO) (70) 
Some comments on the approximations used in this calculation are in order at this point. 

1. We have approximated the vacuum wave function by equation (j36p . This is adequate 
for the linear potential, but there can be corrections due to the higher terms (in powers 
of momenta and powers of J) in the exponent. We hope to address the issue of higher 
order terms elsewhere. 
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2. Wick contractions between W{C) and OjsOa have been neglected. These can lead 
to structures different from W{C) on the right hand side and hence to nondiagonal 
pieces for the first term in (j70p . These are small if the linear potential dominates. The 
lattice data suggest that this is the explained earlier. Further, since W{C) 
has no color charges, these contractions will be suppressed by powers of 

3. We need large M limit to avoid virtual and real processes with additional numbers of 
ip and (p particles. Such processes could also modify the Yang- Mills vacuum, if we do 
not take the large M limit. (This is a quenched approximation for these particles.) 

Consider now the case of the gluelump. To get the wave function of the gluelump, we 
can take equation (|70p with Oq, = 1 and apply on this equation the operator J dJ""q"', where 
7]°- = M^°-^il:^. This has the effect of replacing -0^ by BJ'^M'^"-^ . Further the action of cj)^ on 
the vacuum wave function for the heavy scalars gives 

(f) = Mt-^x- Thus we find 

EE dr{x)W^\x,y)x\y)^Q (71) 
The operator W"''^{x,y) may be written as 

ab 



W''\x,y) 



X 



Pexp / dKK 



—1 



y 



(72) 



It is an open Wilson line for the holomorphic component of the connection dKK ^. (The 
result ()72p follows from the fact that we may think of {A, A) as a complex gauge trans- 
form (by the matrix M^) of {-dHH-^,0), i.e., {A, A) = {-dHH-^,Qf\ and because 
W^\x,y,A3) = g'^{x)W^\x,y)g{y) for A3 = g-^Ag + g^^dg.) 

As in the case of W , W depends on the curve chosen from y to x. For a short separation 
\x — y\, we can write 

W{x,y) = K{x,y)K~\y,y) = U{x,y) (73) 

so that the functional dJ"'W°'^x^ agrees with the gluelump state used in section 4. For a 
longer separation, there are differences. Using the composition property of path-ordered 
exponentials, we find, for a straight-line path, 

W{y + 2e, y) = K{y + 2e,y + e) K-\y + e,y + e) K{y + e, y) K~\y, y) 
^ K{y + 2e, y) K~\y, y) + ee d{dKK~^) + ■■■ 
= U{y + 2e,y) + eed{dKK^^) + --- (74) 
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Thus the Wilson hne W wih differ from U{x,y) = K{x,y)K~^{y,y) by terms involving 
d{dKK~^) and higher derivatives of it when the series is continued. Such terms evidently 
correspond to deformations of the path. 

The gluelump state can be defined, at finite separations, by either dJ"'{x)U"'^{x,y)x^{y) 
or dJ"'{x)W°'^{x,y)x^{y)- Which is more advantageous will depend on which leads to a 
simpler equation within the approximations we are using. The emergence of equation (j70p 
suggests that dJ"-{x)W"'^{x,y)x^{y) is a simpler functional to use, and so we redefine the 
gluelump as 

^G= I d^xcPy f{x,y) dr{x) W''\x,y) x\y) (75) 



Using ([7T]) . we then find 

n j fdJWx ^-0 = [T-C, j dJr]] j /V^^VF</.^^o + (2M + a\x - y\) j fdJWx (76) 

The operator dJrj is a local operator. This must be defined, in a regularized way, by an 
expression of the form Og where f{x,y) is replaced by a{x,y,\); then we must take the 
small A limit. In this case, we can use the result ([56]) : the Coulomb potential is zero and 
we can write 



V2 , 

-M + m ^ ) 

2m 



U^''{u,v)7]\v)^o +••• (77) 



where U"'^{u,v) is used because of the point-splitting needed to define dJr]. The omitted 
terms in ([77]) correspond to terms with either p or p acting on J dJr] and the other left free 
so that it can act on / f^^Wfp^ in ([76]); i.e., they are of the form 

=2 

' Urs{u,v) I [priu), I Bjrj] f7lHPs{v),W]x^ + [Ps{v), I Bjrj] f^Hpriu),W]cl)^ \ 



^ / Urs{u,v) \ [priu), I BJt]] fij^[p,{v),W]<P^ + \ps{v), I BJtj] f ^J^ [pr (u) , W]cP^ \ *o 



e2 



2 



(78) 



One of the terms on the right hand side, namely, e^/2 / Ilrs{u, v)[pr{u), f 9J"] \ps{v), VF"^]x^'' 
is the type of term which led to the Coulomb potential in the calculation of section 4. So one 
might wonder whether this can generate a Coulomb term in addition to the linear potential. 
This is not the case. In fact, the Coulomb potential is due to the term farmBJ"^5{x — u) in 
the commutator [pr{u) , B J"" {x)] = —i&Dar5{x — u). Notice that in (I78p we have the combi- 
nation [pr('u), / Bjrj\ and, since \pr{u)-, r/"(x)] = farmfl^^{x — u) as well, the relevant term is 
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actually zero. (This too has to be understood with point-splitting, but this particular result 
is basically the same as in the unregulated calculation.) We will neglect the remaining terms 
in ()78p . which lead to new operator structures. These terms lead to operator structures with 
more powers of J. The gluelump state we consider can have overlap with states which have 
the same quantum numbers as dJWx- The fact that the Hamiltonian acting on dJWx 
can evolve it into other structures is indicative of this possibility of mixing. However, as is 
well known in quantum mechanics, the effect of such nondiagonal terms in TC comes with 
energy denominators and can be taken to be small if the differences between the energies of 
the lowest gluelump state (which is what we are interested in) and higher states are large 
enough. We expect this to be the case at large enough coupling, since the differences must 
go like m. But ultimately it is to be justified a posteriori^. 

With this understanding of various terms, combining the results in ()76p and ()77p . we 
finally have, as A in ()77p is taken to be very small. 



n 



M + m 



^ X 

2m 



+ a\x - y\ 



f{x,y)dr{x)W-\x,y)x\y)'^o +• 



(79) 



The approximations involved in this equation are given by the comments after equation 
(|70p and by the comments about the omitted terms in (j77l [78|) . 



Equation (j79p shows that the eigenvalue equation we need for the gluelump states is 

V ^ 



M + m 



2m 



+ a\x - y\ 



f{x,y) = E f{x,y) 



(80) 



By our general arguments, the string-breaking energy is given by = 2{Eq—M) where Eq is 
the ground state energy eigenvalue of the equation ([80]) . We now turn to the determination 
of this value. 



6 A variational estimate for string breaking 

Since we are interested in the ground state eigenvalue, the simplest estimate is given by a 
variational calculation. Removing the center of mass motion (which is zero as M ^ oo) the 

^The transition amplitude from the unbroken string state to the gluelumps is seen to be very smaU from 
lattice data, even for SU{2) [121 113j . This is another indication that the off-diagonal elements mentioned 
above might be small. 
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energy functional corresponding to (jSOp is 



E = M + m+ I d^x 



Ml 

2m 



(81) 



We take a variational ansatz of the form / = exp(— where both (3 and ^ are regarded 
as variational parameters. The minimum of ^) with respect to (3 is found to be at 

m/3 



'2(M-3)/M(2ma)r(3//i)V 



A formula for the energy eigenvalue E{P^,iJ,) as a function of /i can now be obtained. 



(82) 



M + m 



2mT{2/fi) \ 



-1/3 



X 8mo-r(3//u) + S^/^/i^ 



'2(A*-3)^2mcjr(3//i)\ 



(83) 



The string tension for a representation R is given by ([2]); for the present calculation, we need 
the adjoint string tension, a a = 7rm^. Using this value of the string tension, the formula 
for E{P^:,iJ,) can be minimized numerically. The minimum is obtained at /x* = 1.752. This 
gives Eq — M = E{P^,iJ,:t:) — M = 3.958m, corresponding to a string-breaking energy of 
Kcaz = 2(S(/3,,/i,) - M) = 7.916 m. 

To check the robustness of the variational result, we have also tried to estimate the 
string breaking energy using a more general variational ansatz of the form 



/ 



(84) 



(1 + \x\y 

Such an ansatz is suggested by the large x solution of the two-dimensional Schrodinger 
equation with a linear potential, whose solutions fall off, at large | 

3/2" 



Ai {{2may/^\x\) exp [-§ {{2may/^\x\) 



„|3/4 



(85) 



where Ai is the Airy function. However, a numerical minimization of i?(/3, /x, z/), with respect 
to the three variational parameters /?, fx, u, shows that the lowest estimate for E is provided 
by = 0, i.e., by the exponential ansatz reported above. 
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We are now in a position to compare this with lattice estimates. The latest data on 
this is for SU{2), given by Kratochvila and Forcrand [13]. The string-breaking energy is 
given by V^^iat ~ 8.68 m. This arises as follows. V^iat is computed in terms of the lattice 
spacmg a as V^iat « 2.063a-^ In the same set-up, the fundamental string tension is given 
by ap = 0.0625a^^. This can be used to write V^iat ~ 2.06 m ^y^^CF/0.0625cA ~ 8.68 m. 
For our analytic calculation, we then have iKiat — Kcaz|/K«at ~ 8.76%. 

It is interesting also to make a similar estimate for the O"*""*" glueball mass, at least as a 
further check on the approximations used. If we use our picture of two constituents of mass 
m each, interacting via the linear potential, the eigenvalue equation is 



2m ^ + a\x -y\ 

2m 2m 



f{x,y)=E f{x,y) (86) 



The only difference, compared to what we have done for the gluelump, is that the reduced 
mass should now be |m. This leads to 

Mo++ = 5.73 m (87) 

The difference with the lattice estimate of 5.17 m [9] is 10.83%. In this case, of course, 
since both constituents have mass m <^ M, our 1/m-expansion and the nonrelativistic 
approximation are less reliable. 



7 Discussion 



A Hamiltonian approach to Yang-Mills theory in 2+1 dimensions was introduced several 
years ago by two of the authors O [H [5]. In this paper we have extended the formalism to 
include matter fields, specifically scalar fields in the adjoint representation. First of all, this 
may be viewed as a first step towards analyzing supersymmetric theories. Secondly, it may 
be used for analyzing screening and string breaking. On general group-theoretic grounds 
we should expect screening of Wilson loops in the adjoint and other representations of zero 
A^-ality. When the length of a string in one of the screenable representations is increased, 
it breaks at some point with the formation of two color singlet states, each of which is 
referred to as a gluelump. We obtained a Schrodinger equation for a gluelump state, the 
variational solution of which then gave an analytic estimate of the string-breaking energy. 
This is within 8.8% of the latest lattice estimates. The emerging picture is that of a bound 
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state between a constituent "gluon" 9J, of mass m, and a heavy adjoint scalar of mass 
M, interacting via a linear potential V ~ (Jj^r. Some approximations are necessary for 
obtaining the Schrodinger equation, since most of the states under discussion are unstable; 
the approximations we have used are given in section 5. It should also be emphasized that 
we are calculating the energy at which the string breaks, the transition amplitude for the 
breaking is not calculated. (These are two different questions. The transition amplitude, 
as is well known, will be suppressed at large N.) 

Specifically for the (2 + l)-dimensional theory, the string-breaking energy has not been 
analyzed extensively on the lattice except for the group SU{2). We consider this to be 
an important question for which more lattice data are needed, for higher groups, and, in 
particular, the A^-dependence. 

Also, as explained in the text, string breaking is closely related to the issue of corrections 
to the formula ([2]) for the string tension. This matter will be taken up in a future publication. 
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Bringoltz. A A also wishes to thank the Max Planck Institut fur Gravitationsphysik, where 
part of this work was carried out, for hospitality. This research was supported in part by the 
National Science Foundation grants PHY-0457304 and PHY-0555620 and by PSC-CUNY 
grants. 

APPENDIX 

Once the regulated Green's functions ()13p have been introduced, it is straightforward to 
evaluate the regularized form of the Hamiltonian and its action on various functionals. 
However, it is important to understand and follow the correct order of limits. Here we shall 
comment briefly on this issue. 

We start by recalling some elementary properties of a covariant integral I{p, e) resulting 
from some term in the perturbation expansion of a field theory. Here p denotes external 
momenta and e is the short-distance cut-off. Once the integration over the loop momenta 
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are done, I{p, e) is of the form 

I{p,e) = J2^"' O^ip) + d{p) + ^e^O^ (Al) 

a p 

The first set of terms involves negative powers of e and represent potentially divergent 
terms as e becomes very small. (We include a = 0, which may be taken as a logarithmic 
divergence.) 0{p) represents finite terms and the last set of terms vanish as e ^ 0. At finite, 
but small, e, the last set of terms are suppressed by powers of p^e; they are negligible for small 
external momenta. While the divergent and finite terms are independent of the regulator 
(apart from the usual ambiguity of the subtraction points), the O'^ terms are sensitive to the 
regulator. Of course, as e ^ 0, this does not matter, but notice that for external momenta 
of the order of e~^/^, these terms can contribute and constitute regularization ambiguities. 
The correct procedure is obviously that the result ()Aip is to be used only for processes with 
external momenta small compared to e~^/^. 

A simple case which highlights this point is the following. A regulator does not have 
to respect symmetries, say, Lorentz invariance; the standard lattice regularization is an 
example. In a Lorentz-invariant theory, this can lead to Lorentz-violating terms of order 
P\/e. If we consistently restrict our calculation to processes with momenta p <C e^^^^, these 
spurious regulator-dependent terms will be negligible. 

In a Hamiltonian formulation, we define the regulated kinetic energy operator T(e), 
given by (jl2p . ()13p . Since it is a functional differential operator, to interpret it properly, we 
have to consider its action on different functionals. Consider a functional ^'(A') with the 
average separation of points between fields in the functional being \/~\! . When we act on 
this with T(e), there are three types of terms possible. 

1. Terms which diverge as e — > 0. If there are such terms, we must introduce a subtraction 
procedure to make T[e) well-defined and eliminate these terms. (The redefinition of 
r(e) with the addition of Q in ()23p may be viewed as an example of this.) 

2. Terms which are finite as e — > 0. These are physical effects. 

3. Terms which are of order e. These can be neglected. However, if we take A' ~ e, 
or A' ^ at finite e, these terms can give finite results, just as the O'^ terms in 
()Aip . Exactly as in the previous case, the correct procedure is to realize that the 
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T(e) can only be used for functionals with A' ^ e; otherwise there will be spurious 
regulator-dependent terms. 



Keeping this in mind, all the calculations have to be done with e ^ first. A natural 
question is then: How do we define a "local" operator like Vym ~ / Note that this 
is needed for the Hamiltonian, and further we will need the commutator of T with this 
operator. From what was said before, Vyai must be defined in a regulated way with an 
average point-separation of \/A, with A ^ e. Then the action of T{e) on this is correctly 
obtained. To make the operator local, we then take A — > 0. In other words, we must take e, A 
to zero, keeping A S> e. This is the correct procedure and it will yield unambiguous results. 
This procedure is even more crucial in a Hamiltonian framework; otherwise, spurious terms 
may appear and it may not even be obvious that they are spurious, since Lorentz symmetry 
is not manifest. 

This procedure was used in 01 [5]. For the purposes of this paper, we note that an 
expression for A^a was obtained in [4J as 

Ara{w, Z) = . ^ . [6ra " {K{w)K-\u, w)K{u)K-\z, u))^ e'^/^j + Ka{w, z) (A2) 
-K^Z — W) /ra J 

A^Jw,z) = -e-''/^\-2d,^—dJK(w)K-^(u,tv)K(u)K-\z,u)) 

IT Z — W ' 

+^—djK(w)K-\u,w)K(u)d,K-\z,u)) 
z — w 

+7 '^——K{w)K-\u,w)K{u)d,{K'\z,u)K{z))K-\z)]^^ 

{Z — W){Z — W) I ira 

+0(e2) (A3) 

where u = \{z^w\ n = ^(z-|-?D) and a = (z — — ?D)/e. If we take to be very small 
compared to the average separation of fields in functionals on which T acts, Ara(it', z) does 
not contribute and we can use the simplified form for Af.a(i?, namely, expression ()A2p 
with Ara(u?, z) set to zero. In fact we can simplify the first term even further, approximating 

kra(w^ Z) as 

Kaiw.z) = ^ [5ra-{K{w)K-\z,w))^^e-''/^] 

■K[Z — W) ™ 

-|- (terms of higher order in e or (z — w), (z — w)) (A4) 

This is correct and adequate, but it must be kept in mind that T so defined acts on operators 
with a point-splitting separation -v/A S> \/e, and in taking the small A-limit we must preserve 
this inequality. 
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